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Abstract 

In this paper we consider a multi-dimensional wave equation with dynamic 
boundary conditions, related to the Kelvin- Voigt damping. Global existence 
and asymptotic stability of solutions starting in a stable set are proved. Blow 
up for solutions of the problem with linear dynamic boundary conditions with 
initial data in the unstable set is also obtained. 
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1 Introduction 



In this paper we consider the following semilinear damped wave equation with 
dynamic boundary conditions: 



Utt — Au — aAut — \u\P 
u{x, t) — 0, 

uttix,t) - 



du , , adut , , , ,,n 9 
— {x,t) + —^{x,t)+rutr~^ 
av ov 



Ut{x,t) 



u{x,0) = ua{x), ut{x,0) = ui{x) 



X eft, t > 
X eVo, t>0 

xeTi, t>0 

X eil , 



(1.1) 



where u = u{x,t) , t > 0, x € , A denotes the Laplacian operator with respect to 
the X variable, is a regular and bounded domain of , {N > 1), dfl = Fq U Fi, 

d 

mes(Fo) > 0, Fq H Fi ~ and — denotes the unit outer normal derivative, 

m > 2 , a , a and r are positive constants, p > 2 and uq , ui are given functions. 

From the mathematical point of view, these problems do not neglect acceleration 
terms on the boundary. Such type of boundary conditions are usually called dynamic 
boundary conditions. They are not only important from the theoretical point of view 
but also arise in several physical applications. For instance in one space dimension, 
the problem (jl.ip can modclize the dynamic evolution of a viscoelastic rod that is 
fixed at one end and has a tip mass attached to its free end. The dynamic boundary 
conditions represents the Newton's law for the attached mass, (see [31 [TJ [3] for 
more details). In the two dimension space, as showed in |26j and in the references 
therein, these boundary conditions arise when we consider the transverse motion of 
a flexible membrane f2 whose boundary may be affected by the vibrations only in a 
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region. Also some dynamic boundary conditions as in problem (jl.ll) appear when 
we assume that fl is an exterior domain of R'^ in which homogeneous fluid is at rest 
except for sound waves. Each point of the boundary is subjected to small normal 
displacements into the obstacle (see p for more details). This type of dynamic 
boundary conditions are known as acoustic boundary conditions. 

In one space dimension, in the case where r 7^ and m = 2, Pellicer and Sola- 
Morales [24j considered the one dimensional problem as an alternative model for 
the classical spring-mass damper system, and by using the dominant eigenvalues 
method, they proved that for small values of the parameter a the partial differential 
equations in the problem (|l.ip has the classical second order differential equation 

miu"{t) + diu'{t) + kiu{t) ^ 0, 

as a limit where the parameter mi , di and ki are determined from the values of 
the spring-mass damper system. Thus, the asymptotic stability of the model has 
been determined as a consequence of this limit. But they did not obtain any rate 
of convergence. This result was followed by recent works [531 [5S]. In particular in 
[25j , the authors considered a one dimensional nonlocal nonlinear strongly damped 
wave equation with dynamical boundary conditions. In other word, they looked to 
the following problem: 

utt - Uocx - autxx + ef (u(l, t), ^ = 0, 

"(0,t)-0, (1.2) 
utf (1, t) = -e [ux + autx + rut] (1, t) - ef (u{l, t), ^ , 

with X E (0, 1), t > 0, r, a > and e > 0. The above system models a spring-mass- 
damper system, where the term ef ^u(l, t), ^ represents a control acceleration 

at a; = 1. By using the invariant manifold theory, the authors proved that for small 
values of the parameter e, the solutions of (11.21) are attracted to a two dimensional 
invariant manifold. See [25], for further details. 

We recall that the presence of the strong damping term —Aut in the problem 
(|l.ip makes the problem different from that considered in [10] and widely studied 
in the literature [301 EH [2H1 El El] for instance. For this reason less results were 
known for the wave equation with a strong damping and many problems remained 
unsolved. Especially the blow-up of solutions in the presence of a strong damping 
and a nonlinear boundary damping at the same time is still an open problem. In 
[llj . the present authors showed that the solution of (|1.1|) is unbounded and grows 
up exponentially when time goes to infinity if the initial data are large enough. 

Recently, Gazzola and Squassina 9 studied the global solution and the finite 
time blow-up for a damped semilinear wave equations with Dirichlet boundary 
conditions by a careful study of the stationary solutions and their stability using 
the Nehari manifold and a mountain pass energy level of the initial condition. 

The main difficulty of the problem considered is related to the non ordinary 
boundary conditions defined on Fi. Very little attention has been paid to this 
type of boundary conditions. We mention only a few particular results in the one 
dimensional space and for a linear damping i.e. (m = 2) [T3l [24l iTllTS]. 

A related problem to (|1.1|) is the following: 

Uu-Au + g{ut) = / mnx{0,T) 

Ou 

— +K(u)utt + h(ut) = 0, ondnx(0,T) 
av 

u{x, 0) — uq{x) in 57 

ut{x,0) — ui{x) in Q, 
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where the boundary term h(ut) — \ut\Put arises when one studies flows of gas 
in a channel with porous walls. The term uu on the boundary appears from the 
internal forces, and the nonlinearity K{u)utt on the boundary represents the internal 
forces when the density of the medium depends on the displacement. This problem 
has been studied in [71 [8]. By using the Faedo-Galerkin approximations and a 
compactness argument, they proved the global existence and the exponential decay 
of the solution of the problem. 

We recall some results related to the interaction of an elastic medium with rigid 
mass. By using the classical semigroup theory, Littman and Markus established 
a uniqueness result for a particular Euler-Bernoulli beam rigid body structure. They 
also proved the asymptotic stability of the structure by using the feedback boundary 
damping. In |19) the authors considered the Euler-Bernoulli beam equation which 
describes the dynamics of clamped elastic beam in which one segment of the beam 
is made with viscoelastic material and the other of elastic material. By combining 
the frequency domain method with the multiplier technique, they proved the expo- 
nential decay for the transversal motion but not for the longitudinal motion of the 
model, when the Kelvin- Voigt damping is distributed only on a subinterval of the 
domain. In relation with this point, see also the work by Chen et al. [5] concerning 
the Euler-Bernoulli beam equation with the global or local Kelvin- Voigt damping. 
Also models of vibrating strings with local viscoelasticity and Boltzmann damping, 
instead of the Kelvin- Voigt one, were considered in [20] and an exponential energy 
decay rate was established. Recently, Grobbelaar-Van Dalsen [T^] considered an ex- 
tensible thermo-elastic beam which is hanged at one end with rigid body attached 
to its free end, i.e. one dimensional hybrid thermoelastic structure, and showed 
that the method used in [21] is still valid to establish an uniform stabilization of the 
system. Concerning the controllability of the hybrid system we refer to the work by 
Castro and Zuazua [4] , in which they considered flexible beams connected by point 
mass and the model takes account of the rotational inertia. 

In this paper we consider the problem ()1.1|) where we have set for the sake of 
simplicity a — \. We will show that if the initial data are in the "stable set", 
the solution continues to live there forever. In addition, we will prove that the 
presence of the strong damping forces the solution to go to zero uniformly and 
with an exponential decay rate. To obtain our results we combine the potential 
well method with the energy method. We will also proved that in the absence of 
a nonlinear boundary damping, in the case where to = 2, the solution blows up in 
finite time. 

This paper is organized as follows. In the section 2, after having stated the 
local existence and uniqueness theorem obtained by the authors in we will 
prove that if the initial data are in the stable manifold, the solution continues to 
live there and so we will prove the global existence and the exponential decay of 
the solution. Let us mention, that despite the Lyapunov functional method is a 
well-know tool to prove exponential decay, at our knowledge, such results were 
not already studied on semilinear damped wave equation with dynamic boundary 
conditions. In section 3, we prove the blow up result of the problem (|l.ip . in the 
case of a linear boundary damping, in spite of the presence of the strong damping 
term Aut. The technique we use follows closely the method used in which is 
based on the concavity argument due to Levine [TB]. Let us mention, also at this 
point that despite the concavity method due to Levine is a well-know tool to prove 
blow up, at our knowledge, such result are new for semilinear damped wave equation 
with dynamic boundary conditions. 
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2 Asymptotic stability 



In this section, we will first recall the local existence and the uniqueness result of 
the solution of the problem (|l.ip porved in [TT] . Then we state and prove the global 
existence and exponential decay of the solution of problem In order to do 

this, a suitable choice of the Lyapunov functional will be made. 
Let us firstly define the critical Sobolev exponent for the trace functional space by: 

+00, if iV= 1,2 



Theorem 2.1 [TT] Let 2 < p < q and max ( 2, 



q 



q+l~p 



< m < q. 



Then given uq G (fl) and ui G L^{fl). there exists T > and a unique solution 
u of the problem lil.l]) on (0, T) such that 



Ut 



G c([o, nnl^^in)) nc' ([o, T],L\n) 

G (0, T; i/i^ nL" ((0, T) X r 



We proved this theorem by using the Faedo-Galerkin approximations and the well- 
known contraction mapping theorem. 

Let us recall that the solution u of (|l.ip belongs to the space: 



Yt = 



Vt G L2(0,r;i7i^(f^))nL™((0,T) x Fi 



endowed with the norm: 



\\u\\y = max 

^ 0<t<T 



\\vt\\l + \\Vv\\ 



L"M(0,T)xri 



\\7vt{s)\\l ds 



Definition 2.1 Let 2 < p < q, max ( 2, ) < m < mq G Hi (fl) and 

V q+^-pJ ° 

Ui G L^(fl). We denote u the solution of We define: 

Tmax — sup|r > , M = u{t) exists on [0, T] | 

Since the solution u £ Yt (the solution is "enough regular"), let us recall that if 
Tmax < 00, then 

lim II Vu||2 + llutlb = +00 . 

If Tmax < oo, we Say that the solution of il.l]) blows up and that Tmax is the blow 
up time. 

If Tmax — 00, we say that the solution of lil.l]) is global. 

In order to study the blow up phenomenon or the global existence of the solution 
of (|1.1|) . we define the following functions: 

m = /Kt)) = ||V^||^-||«||^, (2.2) 

J(t) = J(^^(i)) = ^||Vu||^-i||^||^, (2.3) 
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and 



E{u{t)) = E{t) = J{t) + ]^\\ut\\l + \\\ut\\lr,- (2.4) 



In [llj we obtained the following inequality on the energy 
dE{t) 



dt 



= -a\\^ut\\i-r\\ut\\ll^^<0, Vt>0. (2.5) 



Thus the function E is decreasing along the trajectories. As in [22], the potential 
well depth is defined as: 

d= inf maxJ(Au). (2.6) 

ueHl jn)\{o} A>o 

We can now define the so called "Nehari manifold" as follows: 

M ^ {u e Hl^{n)\{Q}- m^o). 

N separates the two unbounded sets: 

U+ = {ue H^.in); I{t) > 0} U {0} and Af- = {u e H^^{ny,I{t) < O} . 
The stable set W and unstable set U are defined respectively as: 

W = {ue H^^{n);J{t) <d}nAf+ and U = {u e H^^{n);J{t) <d}nAf'. 
It is readily seen that the potential depth d is also characterized by 

d = niin J (u) . 

As it was remarked by Gazzola and Squassina in [5] , this alternative characterization 
of d shows that 



/3 = dist(0,7V) = min ||Vu||2 = W > . (2.7) 

u<^M V p — 2 

In the Lemma [^m we would like to prove the invariance of the set J\f^: if the initial 
data uq is in the set f\f^ and if the initial energy E{0) is not large (we will precise 
exactly how large may be the initial energy), then u(t) stays in A/"^ forever. 

For this purpose, as in P,'^, we denote by C* the best constant in the Poincare- 
Sobolev embedding H^^^{n) ^ LP{n) defined by: 

=M{\\Vu\\2:ueH^^in),\\u\\p^l} . (2.8) 
Let us denote the Sobolev critical exponent: 

_ / ^,ifAr>3 
P= < A-2 ' 

[ +00, if TV = 1,2 

Let us remark (as in [9^1 '30|) that if p < p the embedding is compact and the infimum 
in ()2.8p (as well as in (|2.6p ) is attained. In such case (see, e.g. [22l Section 3]), any 
mountain pass solution of the stationary problem is a minimizer for (|2.8p and C* is 
related to its energy: 

d^P-^C-'^/^"-'^ . (2.9) 
2p 

Let us remark also that in the theorem l2.11 we have supposed that p < q where q is 
defined by (|2.1I) . As q < p, we may use the above characterisation of the potential 
well depth d. 
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Remark 2.1 [351 H] For every solution of (11.11) . given by Theorem 12.11 only one 
of the following assumption holds: 

i) if there exists some to > such that u{to) £ W and E{to) < d, then Vi > 
to , u{t) e W and E{t) < d. 

ii) if there exists some to > such that u(to) G and -^(io) < then Vi > 
to , u(t) e and E{t) < d. 

ui) Vt > , £:(t) > d . 

We can now proceed in the global existence result investigation. For this sake, 
let us state two lemmas. 

Lemma 2.1 Assume 2 < p < q and max ( 2, — — ) < m < q. Let uq G A/"^ 

and ui G L^(r2). Moreover, assume that E{0) < d. Then u{t,.) G Af^ for each 

tG [o,r). 

Remark 2.2 Let us remark, that if there exists to G [0, T) such that 

S(to) < d 

the same result stays true. It is the reason why we choose to = 0. 

Moreover , one can easily see that, from (|2.9p . the condition E{0) < d is equiv- 
alent to the inequality: 

p-2 

Cr(^i?(0)) ' <1 (2.10) 
This last inequality will be used in the remaining proofs. 

Proof of Lemma I2.lt Since I{uo) > 0, then by continuity, there exists T^, < T 
such that I{u{t, .)) > 0, for all t G [0,T*). Since we have the relation: 

P-2„^^„2 , 1 



we easily obtain 



Hence we have: 



^W = ^l|v^||^ + -/(t) 

2p p 



j(t) > ^-^11 v?i||^, vtG[o,r,) 

2p 



p-2 



From and (P^i)) . we obviously have Vt G [0,r*), J(t) < E{t). Thus we obtain: 



\yu\\i < -^E{t) 

p-2 



Since £^ is a decreasing function of t, we finally have: 

_2p_ 

p-2 



Vii||2 < -^£;(0), VtG [0,T,) . (2.11) 



By definition of C*, we have: 



\\u\\i<c:\\vu\\i^c:\\^u\\i-^\\^u\\i 



Using the inequality (|2.1ip . we deduce: 

Ihll? < (^i?(0)) ' liv^ii 



2 
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Now exploiting the inequality on the initial condition (j2.10p we obtain: 

ll«ll^<l|Vu||^ . 

Hence ||Vu||^ - \\u\\p > 0, g [0,T,), this shows that u{t, .) G 7V+, e [0,r,). 
By repeating this procedure, is extended to T. _ 

Lemma 2.2 Assume 2 < p < q and max ( 2, — — ) < m < q. Let uq € A/""*" 



q+l-p, 

andui G L^(r2). Moreover, assume that -B(O) < d. Then the solution of the problem 
is global in time. 

Proof of Lemma 12. 2t Since the map 1 1— 7> E{t) is a decreasing function of time t, 
we have: 

EiO) > E{t) = \\\u,\\l + \\\u,\\lr, + ^l^WVuWl + h{t) , 



which gives us: 



Thus, yt G [0,r) , ||Vu||2 + ||ut||2 is uniformly bounded by a constant depending 
only on E(0) and p. Then by definition I2.1[ the solution is global, so T^ax = oo.^ 

We can now state the asymptotic behavior of the solution of (|l.ip . 
Theorem 2.2 Assume 2 < p < q and max ( 2, - — ) < m < q. Let uq G A/"^ 



q + l-p, 

and ui G L^{fl). Moreover, assume that E{0) < d. Then there exist two positive 
constants C and £, independent of t such that: 

< E{t) < de'^\ V< > 0. 

Remark 2.3 Let us remark that these inequalities imply that there exist positive 
constants K and ^ independent of t such that: 

\\Vu{t)\\l + \\utit)\\l < Ke-^\ yt > 0. 

Thus , this result improves the decay rate of Gazzola and Squassina [9l Theorem 
3.8], in which the authors showed only the polynomial decay. Here we show that 
we can always find initial data satisfying uq G Af^ and ui G L^{fl) which verify 
the inequality (|2.10p . such that the solution can decay faster than 1/t, in fact with 
an exponential rate, even in the case m > 2. Also, the same situation happens in 
absence of strong damping {a — 0) and m — 2. 

Proof of Theorem 12. 2t Since we have proved that Vi > , u(t) G A/"^, we already 
have: 

0<E{t) Vi>0 . 

The proof of the other inequality relies on the construction of a Lyapunov functional 
by performing a suitable modification of the energy. To this end, for e > 0, to be 
chosen later, we define for u G Af~^, 



Vt > , L{t) = E{t) + e 



/ utudx + e / uutda + ^\\\/u\\l . (2.12) 
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Let us see that we have: 



> , \L{t) - E{t)\ 



utudx + e I uutdcr + —\\\/u\\2 



Since we have proved in Lemina [2T1 and Lenima l2.2l that Vt > , I{t) > and ||Vit||2- 
|jMt||2 is uniformly bounded by a constant depending only on E{0) and p, using 
Young inequalities on the two integral terms and then Poincare's inequality, there 
exists a constant C > such that: 



utudx -\- e / uutdfT + — II Vu||2 



< CeE{t) 



Consequently, from the above two inequalities, we have 

> , (1 - Ce)E{t) < L{t) < (1 + Ce)E{t) . 

It is clear that for e sufficiently small, we can find two positive constants j3i and /32 
such that 

Vi > , PiE{t) < L{t) < (32E{t) . (2.13) 

By taking the time derivative of the function L defined above in equation (|2.12l) . 
using problem (jl.ip . and performing several integration by parts, we get: 



^ = -a\\yu,\\l~r\\m\\Z.r,+e\\u^\\l~e\\yu\\l 
+e\\u\\l + £\\ut\\lT^-er I \ut\'^-'^ utuda . 

Now, we estimate the last term in the right hand side of (|2.14p as follows. 
By using Young's inequality, we obtain, for any (5 > 



(2.14) 



\ut\ 



im-2 



Utuda 



< " 

m 



m — 1 



rjj.ri 



m 



(2.15) 



The trace inequality implies that: 

hl|™r, <C||V«r ■ 

where C here and in the sequel denotes a generic positive constant which my change 
from line into line. Since the inequality (|2.1ip holds, we have 



Inserting the two inequalities (|2.15p and (j2.16p in (|2.14l) . we have: 

m — 1 



(2.16) 



'^^^^^ < -a||V.J^ 



dt 



r e 



+e||ut||2+e||ut||2^P^ 



(2.17) 



r8 



-C 



2pEm\ 

p-2 J 



CP 



2p 



* Vb-2) 



E{0) 



V 



From p.lOp . we have 



iVulli 



<o 



CP 



2p 
(p-2) 



£;(o) 



1 < 



8 



Now, let us choose 5 large enough such that: 
P 



p-2 



Once d is fixed, we choose e small enough such that: 



From (|2.17p . we may find 77 > 0, which depends only on 6, such that: 

< -aW^utWl + e\\ut\\l + e\\ut\\lr^ - svW'^uWl 

Consequently, using the definition of the energy (|2.4p . for any positive constant M, 
we obtain: 

^ < -MsE{t) + e[l + ^y\u,\\l^a\\S7u,\\l 

+ (^+s)ll-*lllr.+^(f -^)l|Vu||^ . (2.18) 
By using the Poincare inequality and the trace inequality 

WutWl < cw^utWl 

choosing again e small enough and M < 2ri, from p.lSp . we have: 

dL{t) 



dt 



< -MeE(t), yt > 0. 



On the other hand, by virtue of p.lSp . setting ^ = —Me//32, the last inequality 
becomes: 

^<~m), yt>0 . (2.19) 

Integrating the previous differential inequality (|2.19p between and t gives the 
following estimate for the function L: 

L{t) < Ce-«* , yt>0 . 
Consequently, by using (j2.13l) once again, we conclude 

E{t) < Ce~«* , Vi > . 
This completes the proof. ^ 

Remark 2.4 In [11], we have proved the following result: 

Theorem 2.3 Assume 2 < p < q and m < p. Let uq £ i7p^(ri) and ui S L^(Sl). 
Suppose that 

-p/p-2 



E{0) < d and \\Wuq\\^ > C, 
Then the solution of problem 11. 1]) growths exponentially in the 



norm. 
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The present result on the asymptotic stability completes the above result on the 
exponential growth since when uq e A/''*', we have: ||Vuo||2 ^ C** ^. 

Indeed, since d is the mountain pass level of the function J, we have J(wo) < d. 
This writes: 

P-^\\Vuo\\l + lliO)<d 
zp p 

Since uq e A/"^ , we have: 

n — 9 

Vuo||2<rf . 



2p 

Using identity (EH), we get finally \\Vuo\\^ < C7^^^^'^. 



3 Blow up 

In this section we consider the problem (ll.ip in the linear boundary damping case 
(i.e. m — 2) and we show that ii uq E U and E{Q) < d then any solution of (|l.ip 
blows up in finite time. Our result reads as follows: 

Theorem 3.1 Assume 2 < p < q and m — 2. Let u be the solution of (QHP on 
[0,Tmax)- Then Tmax < oo if and only if there exists t S [OjTmax) such that 

u{t) e U and E{t) < d . (3.1) 

Proof of Theorem 13. It Let us firstly suppose that there exists t e [0, Tmax) 
satisfying the conditions p.ip . Without loss of generality, we may assume that 
t = 0. 

We will prove that Tmax < oo by contradiction. We will suppose that the solution 
is global "in time" and we will use the concavity argument of Levine [16l [17] : the 
basic idea is to construct a positive functional 9{t) of the solution and show that 
for some 7 > 0, the function d~'^{t) is a positive concave function of t. Thus it will 
exists T* such that lim 9~'^(t) ~ 0. From the construction of the function 9, this 

t->T* 

will imply that: 

lim ||Vu||2 + ||ut|l2 = +00 . 

t-j-T* 

t<r* 

In order to find such 7, we will verify that: 



d'9-<{t) 



-/9-''-^{t)99' -{l + j)9"{t) <0, Vt>0 . (3.2) 



dt^ 

Thus it suffices to prove that 9{t) satisfies the differential inequality 

99" ^ {1 + j) 9'\t) > , yt>0 . (3.3) 
From the remark [2. 1[ we firstly have: 

Vt e [0, T,nax) , E{t) < d and u{t) G U . 
Hence by (|2.7p . we obtain: 
2dp 



P- 



<\\Vu{t)\\i, V<e [0,T,„ax) . (3.4) 



Assume by contradiction that the solution u is global "in time". Then for any 
T > 0, let us define the functional 9 as follows 



9{t) = \\u{t)\\l + \\u{t)\\lr^+a \\Vu{s)\\jds + r \\u (s) Wl^^ds 

Jo Jo 

+ {T-t)[a\\Vuo\\l + r\\uo\\lr,], Vte[0,T). (3.5) 
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Taking the time derivative of (13.51) we have: 

0'{t) = 2 Utudx + 2 / Utuda + 2a / / WuWutdxds 
Jn JTi Jo Jn 



'2r 



Ufudads 



(3.6) 



Jr 



Thus, as u is the solution of problem (|1.1[) . differentiating equation p. 61) with respect 
to t gives us: 

9"it) = 2 [Wutiml \\Vuit)\\l + \\u\\; + Wutmlr,] . 
Therefore, using the definition of 9 given by p.Sp . we can easily see that: 
P + 2, 



eit)e (t) 



-0 {tf = 20it) [\\u,{t)\\l - \\Vuit)g + \\u\\P + \\u,{t)\\lr^] 



e{t)-{T-t) a\\Vuo\\j + r\\uo\\lr, 



(3.7) 



lk*Wll2 + ll"iWlllri +a / ||Vut(<)||^ds + r / ||ut(i)||lrirfs 



+ (p + 2) 77(i) 
where the function r/ is defined by: 



vit) 



\Ht)\\2 + \Ht)\\2T,+a \\Vuit)\\ids + r \\uit)\\ir^ds 



\\utit)\\i + \\utit)\\l^^ + a \\Vutit)\\ids + r Wutmirjs 



(3.8) 



Utudx - 



Ufudcr + a 



ri 



VuVutdxds + r 



Jn 



Utudads 



JTt 



Our purpose now is to show that the right hand side of the equality (13. 7p is non 
negative. Let us firstly show that 'q{t) > for every t G [0,T]. To do this, we 
estimate all the terms in the third line of (|3.8p making use of Cauchy-Schwarz 
inequality, and compared the results with the terms in the first and second line in 
p.8p . For instance, when we develop the square term in the inequality (|3.8I) . we 
estimate the terms as follows: 



utudx ) < ||'"(i)|l2ll"t(^)ll2 

2/ utudx I Ufuda < ||u(t)||2,ri Il"t(^)ll2 + Il^'t(*)ll2,ri Il"(^)ll2- 
Jn JTi 

Also, the following estimate holds: 

2a / \7u\/utdxds / utudx < a\\ut{t)\\2 / || Vu (s) ||2(is 
Jo Jn Jn Jo 

+a\\uit)\\l f\\^ut{s)\\lds . 
Jo 

By carrying "carefully" all computations based on the same estimates as above, we 
finally obtain 

Vte[0,r],ry(t)>0 . 
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Consequently, the equality p.7|) becomes 

ditW" (t) - {tf > d{t)C{t) , yt e [0, T] 



where 



at) = 2[h,(i)||^-||V^(t)||^ + h||^ + h,(i)|llrJ 
-ip + 2){\\utit)\\l + \\utmlr, 

+a I \\Vut{t)\\lds + r f \\ut{t)\\l r,ds} 



Let us remark that 

at) = -2pE{t) + {p-2)\\Wu{t)\\l-{p + 2)a f \\Wut{s)\\lds 

-{p + 2)r f \\utis)\\lrjs 







From the equality (|2.5p . we have: 

yte[0,T],Eit)+a f \\Vutis)\\lds + r f \\ut{s)\\lr^ds ^ E^O) . (3.9) 



JO 



Thus we can write: 



at) = -2pE{0) + {p-2)\\Vu{t)\\l-{p-2)\\S/u„\\l 



+ {p~2)a I \\Vut{s)\\lds + {p-2)r I \\ut{s)W'ds. 



Therefore, by using p.4p and since £'(0) < d we have: 



at) > 2pid-E{0)) + ip~2)a f \\^ut{s)\\lds + {p-2)r [ \\ut{s)\\l^ds 

Jo Jo 

> ip-2)a f \\Vut{s)\\lds + {p-2)r f \\utis)\\ij.^ds. 
Jo Jo 

Hence, there exists S > such that 

at)>s, yte[to,T] 

for some to > 0. But, since 0{t) is continuous and positive, there exists p > such 
that 

9it)>p, yte[to,T] 

Consequently, 

9it)e"{t) - ^e'{t)^ >ps, yte [to, t] 

Setting 

p-2 

the differential inequality p.3p is verified on [to, T]. This proves that 9{t)~^ reaches 
in finite time, say as f — T*. Since T* is independent of the initial choice of T, 
we may assume that T* < T. This tells us that: 

lim e(t) = +0O . 
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In turn this implies that: 

lim ||Vu||2 + ||ut|l2 = +00 . 

t<T* 

Thus we cannot suppose that the solution of (|l.ip with rn = 2 is gobal " in time" , 
that is Tmax < 00. 

Conversely, let us suppose that Tmax < 00. 
By Holder's and Young's inequalities, we have: 



2l(i ^""*"2.ri+V^I|Vw*ll2C?.s) <^ {a\\ut\\lr,+r\\yut\\l)ds 



This implies 

ft 







V 2 



2 



Since p > 2, we get form (12.31) and 



Int) = I {\\yuit)\\i \\u{tw,) < Iwvumi - ^hmi < Eit) (3.10) 

By the help of (13. 9p . we thus have: 



{MHo)\kT, + M^mhfy (3.11) 



Since Tmax < oo, the solution is not global in time. Thus by the defintion of T„ 
(see defintion l2.ip . we have: 



Im [\\u{t)\\^^^^ + \\Vu{t)\\^] ^+00 (3.12) 

* J- max ^ 

From ([XTU)) and (jXTTI) . we obtain 



lim I(t) = lim E(t) = -oo. (3.13) 

Since the functions /(t) and E(t) are continuous, there exists ti € [0,T!,„aa.) such 
that 

E{ti) < d and /(ti) < . 
Thus the time ti verifies the condition (13. 1|) . ■ 



Remark 3.1 The term f{u) = \u\p^'^u is clearly responsible for the blow up sit- 
uation. It is often called the "blow up term". Consequently when f{u) = 0, or 
f{u) = —\u\^~'^u any solution with arbitrary initial data is global in time and the 
result of Theorem 13.11 holds without condition (|2.10p . 

Remark 3.2 It's early well known ( |l61 fT7) ) that this blow up result appears for 
solutions with large initial data i.e. E{0) < 0. We note here that if £^(0) < 0, then 
the blow up conditions p.ip hold. 



Acknowledgments 

The second author was partially supported by MIRA 2007 project of the Region 
Rhone-Alpes. This author wishes to thank Univ. de Savoie of Chambery for its 
kind hospitality. 



13 



References 



[1] K. T. Andrews, K. L. Kuttler, and M. Shillor. Second order evolution equations 
with dynamic boundary conditions. J. Math. Anal. AppL, 197(3):781-795, 
1996. 

[2] J. T. Beale. Spectral properties of an acoustic boundary condition. Indiana 
Univ. Math. J., 25(9):895-917, 1976. 

[3] B. M. Budak, A. A. Samarskii, and A. N. Tikhonov. A collection of problems 
on mathem,atical physics. Translated by A. R. M. Robson. The Macmillan Co., 
New York, 1964. 

[4] C. Castro and E. Zuazua. Boundary controllability of a hybrid system con- 
sisting in two flexible beams connected by a point mass. SIAM J. Control 
Optimization, 36(5): 1576-1595, 1998. 

[5] S. Chen, K. Liu, and Z. Liu. Spectrum and stability for elastic systems with 
global or local kelvin-voigt damping. SIAM J. Appl. Math., 59(2):651-668, 
1999. 

[6] F. Conrad and O. Morgiil. On the stabilization of a flexible beam with a tip 
mass. SIAM J. Control Optim., 36(6):1962-1986 (electronic), 1998. 

[7] G. G. Doronin and N. A. Larkin. Global solvability for the quasilinear damped 

wave equation with nonlinear second-order boundary conditions. Nonlinear 
Anal, Theory Methods Appl. 8:1119 1134, 2002. 

[8] G. G. Doronin, N. A. Larkin, and A.J. Souza. A hyperbolic problem with 
nonlinear second-order boundary damping. Electron. J. Differ. Equ. 1998, 
paper 28, pages 1-10, 1998. 

[9] F. Gazzola and M. Squassina. Global solutions and finite time blow up for 
damped semilinear wave equations. Ann. I. H. Poincare, 23:185-207, 2006. 

[10] V. Georgiev and G. Todorova. Existence of a solution of the wave equation with 
nonlinear damping and source terms. J. Differential Equations, 109(2) :295-308, 
1994. 

[11] S. Gerbi and B. Said-Houari. Local existence and exponential growth for a semi- 
linear damped wave equation with dynamic boundary conditions. Advances in 
Differential Equations, 13(11-12):1051 1074, 2008. 

[12] M. Grobbelaar-Van Dalsen. Uniform stabilization of a one-dimensional hybrid 
thermo-elastic structure. Math. Methods Appl Sci., 26(14):1223-1240, 2003. 

[13] M. Grobbelaar-Van Dalsen and A. Van Der Merwe. Boundary stabilization 
for the extensible beam with attached load. Math. Models Methods Appl. Sci., 
9(3):379-394, 1999. 

[14] Esquivcl-Avila J. The dynamics of nonlinear wave equation. J. Math. Anal 
Appl, 279:135 150, 2003. 

[15] M. Kirane. Blow-up for some equations with semilinear dynamical boundary 
conditions of parabolic and hyperbolic type. Hokkaido Math. J., 21(2):221-229, 
1992. 

[16] H. A. Lcvinc. Instability and nonexistence of global solutions to nonlinear wave 
equations of the form Pwtt = —Au + T{u). Trans. Amer. Math. Soc, 192:1-21, 
1974. 



14 



[17] H. A. Levine. Some additional remarks on the nonexistence of global solutions 
to nonlinear wave equations. SIAM J. Math. Anal., 5:138-146, 1974. 

[18] W. Littman and L. Markus. Stabilization of a hybrid system of elasticity by 
feedback boundary damping. Ann. Mat. Pura Appl., IV. Ser. , 152:281-330, 
1988. 

[19] K. Liu and Z. Liu. Exponential decay of energy of the Euler-BernouUi beam 
with locally distributed Kelvin- Voigt damping. SIAM J. Control Optimization, 
36(3):1086-1098, 1998. 

[20] K. Liu and Z. Liu. Exponential decay of energy of vibrating strings with local 
viscoelasticity. Z. Angew. Math. Phys., 53(2):265-280, 2002. 

[21] K. Ono. On global existence, asymptotic stability and blowing up of solutions 
for some degenerate nonlinear wave equations of KirchhofF type with a strong 
dissipation. Math. Methods Appl. Sci., 20(2):151-177, 1997. 

[22] L. E. Payne and D. H. Sattingcr. Saddle points and instability of nonlinear 
hyperbohc equations. Israel J. Ma.th., 22(3-4):273 303, 1975. 

[23] M. Pellicer. Large time dynamics of a nonlinear spring-mass-damper model. 
Nonlin. Anal, 69(1):3110-3127, 2008. 

[24] M. Pellicer and J. Sola-Morales. Analysis of a viscoelastic spring-mass model. 
J. Math. Anal. Appl, 294(2):687-698, 2004. 

[25] M. Pellicer and J. Sola-Morales. Spectral analysis and limit behaviours in a 
spring-mass system. Commun. Pure Appl. Anal, 7(3):563-577, 2008. 

[26] G. Ruiz Goldstein. Derivation and physical interpretation of general boundary 
conditions. Adv. Differ. Equ., ll(4):457-480, 2006. 

[27] G. Todorova. The occurence of collapse for quasilinear equations of parabolic 
and hyperbolic type. C. R. Acad Sci. Paris Ser., 326(1):191-196, 1998. 

[28] G. Todorova. Stable and unstable sets for the cauchy problem for a nonlinear 
wave with nonlinear damping and source terms. J. Math. Anal Appl, 239:213- 
226, 1999. 

[29] G. Todorova and E. Vitillaro. Blow-up for nonlinear dissipative wave equations 
in R". J. Math. Anal Appl, 303(l):242-257, 2005. 

[30] E. Vitillaro. Global nonexistence theorems for a class of evolution equations 
with dissipation. Arch. Ration. Mech. Anal, 149(2):155-182, 1999. 



15 



